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ABSTRACT 

We use a recently released SDSS catalog of X-ray emitting AGNs in con- 
junction with the FIRST 20cm radio survey to investigate the black hole fun- 
damental plane relationship between the 1.4GHz radio luminosity (L r ), 0.1-2.4 
keV X-ray luminosity (Lx), and the black hole mass (M), namely, logL r = 
£,RxlogL x + £,RMlogM+const. For this purpose, we have compiled a large sample 
of 725 broad-line AGNs, which consists of 498 radio-loud sources and 227 radio- 
quiet sources. Our results are generally consistent with those in our previous 
work based on a smaller sample of 115 SDSS AGNs. We confirm that radio-loud 
objects have a steeper slope (£_rx) in the radio- X-ray relationship with respect 
to radio-quiet objects, and the dependence of the black hole fundamental plane 
on the black hole mass (£rm) is weak. We also find a tight correlation with a 
similar slope between the soft X-ray luminosity and broad emission line lumi- 
nosity for both radio-loud and radio-quiet AGNs, which implies that their soft 
X-ray emission is unbeamed and probably related to the accretion process. With 
the current larger sample of AGNs, we are able to study the redshift evolution 
of the black hole fundamental plane relation for both radio-loud and radio-quiet 
subsamples. We find that there is no clear evidence of evolution for radio-quiet 
AGNs, while for radio-loud ones there is a weak trend where ^rm decreases as 
the redshift increases. This may be understood in part as due to the observed 
evolution of the radio spectral index as a function of redshift. Finally, we discuss 
the relativistic beaming effect and some other uncertainties related to the black 
hole fundamental plane. We conclude that, although introducing scatters to the 
fundamental plane relation, Doppler boosting alone is not enough to explain the 
observed steeper value of £rx in the radio-loud subsample with respect to the 
radio-quiet ones. Therefore, the significant difference of £rx between radio-loud 
and radio-quiet sources is probably also due to the different physical properties 
of the jets. 
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1. INTRODUCTION 

Astrophysical black holes do not emit light directly but they can be probed by their 
gravitational influence on neighboring matter, which produces observable signatures of black 
hole activity. The key mechanism for a black hole to become active is the accretion process 
(Frank, King & Raine 2002), which is usually accompanied by a relativistic jet. Accretion 
disks and jets can produce photons from radio to X-ray band. The radio emission is usually 
believed to originate from the synchrotron radiation of the jet (Begelman, Blandford & Rees 
1984), while the optical/UV emission mostly comes from the multicolor black body radiation 
emitted from the accretion disk (Shakura & Sunyaev 1973), and the X-ray radiation is usually 
associated with the inner-most region of the accretion disk, where the temperature is the 
highest. In some cases, the contribution of inverse Compton scattering from high energy 
electrons in a disk corona is also needed to account for the observed power-law spectrum 
in the X-ray band (Haardt & Maraschi 1993). If jet production is directly related to the 
accretion process, we would expect a natural correlation between radio and X-ray luminosities 
(Merloni, Heinz & Di Matteo 2003; Heinz & Sunyaev 2003; Falcke et al. 2004). 

The radio-to- X-ray correlation has long been studied in both Galactic black hole (GBH) 
candidates and active galactic nuclei (AGNs). Gallo, Fender & Pooley (2003) found a strong 
correlation between the radio and X-ray emission (L r oc L ^) using the simultaneous X- 
ray and radio observational data of stellar-mass black hole X-ray binaries (XRBs) during 
the low/hard state. In addition, they suggested that when XRBs enter the hard to soft 
transition state, the jet is suppressed and the radio emission decreases. Recently, some 
studies have shown that the substantial scatter exists in such a relationship of GBHs (Gallo 
2006; Xue & Cui 2007; Xue, Wu & Cui 2008). In the case of AGNs, Brinkmann, Yuan & 
Siebert (1997) obtained a remarkable correlation between 2keV X-ray luminosity and 5GHz 
radio luminosity for 324 radio-loud AGNs. Canosa et al. (1999) found a strong correlation 
between soft X-ray and 5GHz radio luminosities for 40 low-power radio galaxies. Brinkmann 
et al. (2000) studied a sample derived from the cross- correlation of the ROSAT All-Sky 
Survey (RASS) catalog and the Very Large Array (VLA) Faint Images of the Radio Sky at 
Twenty-Centimeters (FIRST 20cm) catalog. They found that for 843 AGNs, the X-ray and 
radio luminosities are correlated over two decades in radio luminosity, spanning radio-loud 
and radio-quiet regimes, but radio-quiet quasars seem to follow a different correlation from 
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radio-loud ones. Recently, Panessa et al. (2007) investigated the radio/X-ray luminosity 
correlation for low-luminosity AGNs, including local Seyfert galaxies and low luminosity 
radio galaxies (LLRGs). They found that X-ray and radio luminosities are significantly 
correlated over 8 orders of magnitude, with Seyfert galaxies and LLRGs showing the similar 
slope, which seems different from the previous results. 

Theoretical explanations for the observed radio-to-X-ray relation have been discussed 
in a number of previous works. Fender et al. (2003) found that, at relatively lower accretion 
rate ( M < 7 x 10^ 5 M^^), black hole X-ray binaries would enter a "jet-dominated" state. 
At this stage, the majority of the liberated accretion power is transferred into the jet and 
does not dissipate as X-ray emission in the accretion flow. This was also suggested by Falcke, 
Kording & Markoff (2004), who demonstrated that, below a critical value of the accretion 
rate (< 1 ~ 10%M Edd ), the Spectral Energy Distribution (SED) of a black hole accreting 
source is dominated by the non-thermal radiation from the jet, while for sources with higher 
accretion rates (M < MEdd), the radiation is dominated by the accretion flow. Using a 
coupled disk-jet model, Yuan & Cui (2005) showed that, in GBHs when the X-ray luminosity 
is smaller than a certain value (~< 10~ 5 — lO^L^d), the jet will dominate the radiation and 
the radio-to-X-ray correlation becomes to L r oc L x 23 - When the X-ray luminosity exceeds 
that critical value, the accretion flow would produce most of the X-ray emission and the 
radio-to-X-ray correlation then becomes to L r oc L x 7 - Therefore, the accretion rate decides 
the detailed physical model of the accretion disk, which then leads to different observed 
radio-to-X-ray correlation slopes (Kording, Falcke & Corbel 2006). Although it is still not 
clear whether such a critical accretion rate is of the same order of magnitude for both GBHs 
and AGNs, the similarity of accretion-jet physics in these two systems seems to imply that 
similar results may also be found in AGNs. 

In a comprehensive study, Merloni et al. (2003) examined a sample combining Galactic 
black hole systems and supermassive black hole systems by investigating their compact 
emission in X-ray (2-10 keV) and radio (5 GHz) bands. They found that the radio luminosity 
(L r ) is strongly correlated with both the black hole mass (M) and the X-ray luminosity 
{Lx) (so called the fundamental plane of black hole activity). The relation is: logL r = 
(0.601^) logL x + (0.78t f 9 ) logM+ (7.33l£oy)- Subsequently, Wang et al. (2006) selected 
a uniform sample of broad-line AGNs which was cross-identified from the RASS, Sloan Digital 
Sky Survey (SDSS) and FIRST 20cm radio survey to test the black hole fundamental plane 
relation. Their final sample consisted of 115 broad emission line AGNs including 39 radio- 
quiet AGNs and 76 radio-loud ones. They found that the relationship has a very weak 
dependence on the black hole mass. Moreover, radio-quiet and radio-loud objects have 
different radio-to-X-ray slopes, which is 0.85 for radio-quiet objects and 1.39 for radio-loud 
sources. This differs from the result of Merloni et al. (2003) where the relationship seems 
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to be universal for different types of black hole sources. However, the limited statistics of 
the sample in Wang et al. (2006) motivated us to increase the number of sources in order to 
confirm the results on stronger statistical bases. With a larger sample, we are also able to 
study the possible evolution of the black hole fundamental plane relation as we have enough 
sources in each redshift bins. 

We organize the paper as the following. In §2 we present our sample selection criteria and 
the main properties of our sample. In §3 we show the derived fundamental plane relation and 
investigate the possible evolution of such a relation. In §4 we briefly discuss and summarize 
our results. 



2. THE RADIO AND X-RAY EMITTING BROAD-LINE AGN SAMPLE 



2.1. Sample selection 

Our sample is selected based on the cross-ide ntification of the newly published X-ray 



emitting SPSS AGN catalo g 
gjT fTWhite et al.l [l997h . 



surve\ 



Anderson et al.ll2OO70 and the catalog of the FIRST 20cm radio 



Anderson et al.l (120071 ) employed X-ray data from the ROSAT 



All-Sky Survey (RASS) and both optical imaging and spectroscopic data from the Sloan 
Digital Sky Survey (SDSS). It is worth emphasizing that the RASS and SDSS are extremely 
well matched to each other via a variety of coincidences (e.g., similar survey depth, sensitivity 
and sky coverage, etc.) (Anderson et al. 2003). The RASS/SDSS data from 5740 deg 2 of sky 
spectroscopically covered in SDSS Data Release 5 provide an expanded catalog of 7000 



confi rmed quasars and other AGNs that are probable RASS identifications (lAnderson et al. 
2007h . 



One of the main benefits of this sample is that the SDSS surveyed area is also cov- 
ered by the FIRST 2 0cm radio survey. We cross-correlate the broad-line AGN catalog of 
Anderson et al.l (120071 ) with the FIRST radio-detected sources to build a RASS-SDSS-FIRST 
cross-identified sample of 868 broad-line AGNs. All of these 868 sources have been observed 
and detected at 1.4GHz (FIRST 20cm survey) and in 0.1-2.4 keV energy band (RASS). 
The optical spectra from the SDSS data archive can be used to estimate the central black 
hole mass. Here, we exclude the high z sources (z > 2.171) in our sample for the lack of 
the H/? and Mg II A2798 emission lines in their SDSS optical spectra. Black hole masses 
of the remaining 725 sources are estimated. The radio-loudness (R) is calculated with the 
rest-frame B-band (4400A) and 5GHz flux density according to the definition R = f^GHz/ fs 



1 See VizieR Online Data Catalog, 8071 (Becker et al. 2003) 
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(Kellermann et al. 1989), and the 5GHz flux density is derived from the 1.4GHz flux den- 
sity assuming a spectral index of 0.5 (Kellermann et al. 1989). Radio-loud and rad io-quiet 
sources are separated by R = 10. Please see subsection 2.2 of 1 Wang et al.l (120061 ) for the 
details of data reduction and black hole mass estimation (Kaspi et al. 2000, 2005; McLure 
& Javis 2002; Wu et al. 2004). 

Our final sample comprises 725 entries, of which 498 are radio-loud, 227 are radio-quiet. 
In addition, we discover several radio-loud narrow-line Seyfert 1 galaxies (NLS1) in our study. 
A detailed discussion of these radio-loud NLSls is presented in Appendix A. Throughout 
this paper, we use the cosmologj@ with H = 70 kms _1 Mpc _1 , Ha = 0.7 and Qm = 0.3. 



2.2. Sample properties 

We describe here the main properties of our selected AGN sample. Table [T] gives the 
total 725 sources with the SDSS optical name, redshift, broad-band (0.1-2.4 keV) soft X- 
ray luminosity, rest-frame 1.4GHz radio luminosity, radio- loudness, black hole mass, broad 
emission line luminosity and corresponding flag. In Figure [U histograms of the redshift, 
logarithm of radio-loudness, 0.1-2.4 keV X-ray luminosity, black hole mass (in unit of M ), 
1.4GHz radio luminosity and ratio of X-ray to Eddington luminosity (Le^jj of our sample 
are plotted. The radio-loudness distribution does not show a clear dichotomy between radio- 
loud and radio-quiet AGNs. This is not surprised for samples of FIRST detected quasars. 
The radio-loudness dichotomy was found mostly for optically selected samples, and is usually 
absent for samples of deep radio-detected sources (Brinkmann et al. 2000; White et al. 2000; 
Lacy et al. 2001). 

Our sample is a part of the broad-line AGN catalog provided by Anderson et al. (2007), 
which includes typical broad-line quasars, Seyfert 1 galaxies, low redshift Seyfert 1.5-1.9 
galaxies and some rare galaxies like NLSls. Compared with Wang et al.(2006), the dis- 
tribution ranges of different physical parameters in our sample are relatively larger. The 
significant range in the luminosity at each redshift bin can avoid the strong dependence of 
luminosity on redshift in a flux-limited sample (Avni & Tananbaum 1982). In particular, 



2 We now use this new cosmology to calculate the black hole mass from the H/3 broad 
component's luminosity and FWHM. That is, M bh = 1-464 x 10 5 (y^ 8 ^) ( in™"s ) and 
LogR B LR = (1-324 ± 0.086) + (0.667 ± 0.101)Log ( wii ^ s _ t ) . The latter relation is an updated 
version of the Rblr — Lh/j relation proposed by Wu et al. (2004) using the new cosmology. 

3 Defined as L Edd = 1.26 x 1O 38 (M/Af ) ergs s" 1 
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more convincing results can be obtained with the larger sample. 

However, we must notice the selection biases in our sample. First of all, AGNs with 
radio flux fainter than the FIRST detection limit can not be detected by FIRST. Therefore, 
our sample does not include radio-quiet AGNs with logio (Radio — Loudness) < -0.230 In 
addition, based on the redshift distribution shown in Figure 1, it is clear that most AGNs at 
high redshifts in our sample are radio-loud. Therefore, any conclusion based on the total sam- 
ple is biased toward luminous radio-loud AGNs. We also perform the Kolmogorov-Smirnov 
(K-S) test to evaluate the distribution similarity of physical parameters shown in Figure 1 
between radio-quiet and radio-loud subsamples (such as L r , Lx, M BH , and ^-). We found 
that the probabilities of similarity are all less than 0.05, meaning that the distribution of 
these physical quantities in the two subsamples are essentially different, with relatively larger 
mean value of each parameter for radio-loud AGNs. 

We plot the radio vs. X-ray luminosity (left panel) and the Eddington-luminosity-scaled 
radio vs. scaled X-ray luminosity (right panel) in Figure [2j Objects in different black hole 
mass bins are presented with different symbols and colors to highlight a possible segregation 
in the plot. Clearly, there is no trend that sources in different black hole mass bins are 
parallel to each other, which is consistent with the result in Wang et al. (2006). Again, tight 
correlation between the radio and X-ray luminosities is clear, with radio luminosity spanning 
more than 6 orders of magnitude. 

Figure [3] and Figure 2 are identical except that in Figure 3, different colors and symbols 
are used to denote different radio-loudness bins instead of black hole mass bins. AGNs in 
different radio-lo udness bins seem to distribute in parallel sequences. The result also confirms 



the conclusion of IWang et al.l (120061 ). 



3. CORRELATION ANALYSIS 

In this section, we derive the black hole fundamental plane relation based on the current 
large sample. We also adopt several statistical methods to test its significance. Finally, we 
examine the possible evolution of this relation at different redshifts. 



10 °' 23 is the minimum radio-loudness value of radio-quiet AGNs in our sample 
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3.1. The fundamental plane relation 

In the three-dimensional space (logL r , logLx, logM), black hole systems ar e preferentially 



distri buted on a plane, called the "fundamental plane of black hole activity" (IMerloni et al. 



20031 ) . Based on our sample of AGNs, we fit the data in the form: 

l0g (lO-ergs 8-0 = eRxl ° g (lO-ergss-0 + W ° g (l0^) + ^ (1) 

We also fit the relation between Eddington-luminosity-scaled radio and X-ray luminosi- 
ties for the radio-quiet subsample only. The fitting formula is 

log f J = £ E Rxlog ( 7—^- J + constant. (2) 



\LEdd/ \LEdd/ 

We apply the Ordinary Least Squares(OLS) multivariate regression method (Isobe et 
al.1990) to the total, radio-loud and radio-quiet subsamples, respectively. The OLS bisector 
fitting result for equation (1) with errors at the one-sigma confidence level and the dispersion 
(ovjj are given in Table [2J We also list results of previous works for comparison. Our result 
is generally consistent with Wang et al. (2006), but now with smaller uncertainty for each 
coefficient due to the larger sample. The conclusions are similar: first, the coefficient £rx 
tends to be larger as the radio- loudness increases; second, the black hole mass seems to be 
unimportant in the correlation between radio and soft X-ray luminosities. We should note 
that the black hole fundamental plane coefficients calculated by Merloni et al. (2003) are 
based on a radio-loud and radio-quiet combined sample, with a predominance of radio-quiet 
sources. Our result for the radio-quiet subsample is similar to that of Merloni et al. (2003) 
but the dependence on the black hole mass is weaker in our case. 

The fitting result for equation (2) with errors at the one-sigma confidence level is 

log ( — ^- J = (0.96 ± 0.04) log ( ~~~ J + (-4.82 ± 0.08) . (3) 

\ LEdd / \ LEdd / 

The dependence of L r /Lx on the black hole mass is shown in Figure 4. The overall corre- 
lation is weak, except that a positive correlation is observed for sources (mostly radio-quiet 
ones) with black hole masses smaller than 10 7 ' 5 M Q . Figure 5 shows the edge-on black hole 
fundamental plane relation for radio-quiet sources, with radio-loud AGNs overplotted for 
comparison. We also plot L r /L E dd vs. L x /L E dd in Figure 6, with different symbols represent- 
ing different radio-loudness bins. Parallel sequences can be seen clearly from these figures. 



5 We define the dispersion as the square root of the variance of the differences between the observed radio 
luminosity and that calculated from the fitting relation 
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3.2. Statistical tests for the fundamental plane relation 



3.2.1. Partial correlation tests 



As lBregman et al.l (120051 ) pointed out, the correlation between X-ray and radio luminosi- 
ties may be dominated by the distance effect. Following IWang et al.l (120061). we perfor med 



the partial Kendall r correlation test to examine this effect (lAkritas fc Siebert 1 119961 ). In 
Table El Column (l)-(3) list the variable names of X, Y and Z respectively, where the par- 
tial correlation of X and Y is calculated with the influence of Z variable excluded. Column 
(4) gives the subsample type. Column (5) lists the number of sources in each subsample. 
Column (6)- (8) show results of the partial correlation test, the square root of the calculation 
variance and the probability of the null hypothesis. The null hypothesis is rejected with a 
probability less than the significance level(i.e.,~ 0.05). From Table [3] we can see that the 
partial correlation between L r and Lx is strong even after excluding the distance effect, 
because the P nu ii value is less than 10" 
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3.2.2. The scrambling test 

Besides the partial correlation test performed above, we adopt another method intro- 
duced by Bregman (2005) to evaluate the degree of influence that any distance effect has 
on our sample (i.e., Merloni et al. 2006). We calculate the Pearson correlation coefficient 
(p) between Log ( 10 4o e ^ and ZnxLog ( j^f s pr ) + ^RAfLog ( 10 t\J by randomly 
assigning radio fluxes to objects in our sample. This procedure is performed 10 6 times in 
order to construct the Monte Carlo test. The result is shown in Figure 7. The correlation 
coefficients are adopted from Table 2. For comparison, we overplot the Pearson correlation 
coefficient of the original sample with a vertical line in each panel of Figure 7. 

For radio-quiet objects, our 10 6 realizations of randomized datasets produce only one 
case, where the Pearson correlation coefficient exceeds the value of the real datasets. For 
radio-loud sources, not even one shows a stronger correlation than the real value. This means 
that, for the radio-quiet subsample, the probability that the correlation of the fundamental 
plane is entirely due to the distance effect is about 1 x 10 -6 , and it is even less than 10~ 6 
for the radio-loud subsample and the entire sample. Therefore, we are confident to say that 
the existence of a correlation between the radio luminosity, X-ray luminosity and black hole 
mass cannot be purely the result of distance effects. This result is consistent with the partial 
correlation test performed above. 



- 9- 



3.3. The correlation between soft X-ray and broad emission line luminosities 

In Figure 8, we plot the soft X-ray luminosity versus the H/3 broad emission line lumi- 
nosity (left panel) and the radio luminosity versus the the Hf3 broad emission line luminosity 
(right panel) in our sample. Radio-loud and radio-quiet sources are marked with different 
symbols to highlight the possible dichotomy of the slope in the figure. It is apparent in 
Figure 8 that radio-loud and radio-quiet objects almost have the same slope in the left panel 
(L x vs. L H p). The slope is 1.04 ±0.04 for the radio-quiet subsample and 1.00 ±0.04 for the 
radio-loud subsample. The Pearson correlation coefficients are 0.86 and 0.80, respectively. 
However in the right panel of Figure 8 (L r vs. L H/3 ), the slope of radio-loud sources is 
steeper than that of radio-quiet ones. The slope is 1.68 ± 0.08 for radio-loud sources and 
0.99 ±0.03 for radio-quiet sources, and the Pearson correlation coefficients are 0.64 and 0.88, 
respectively. 

The similar slope of radio-loud and radio-quiet sources in the Lx vs. Lup plot seems 
to indicate that the soft X-ray emission traces well the ionizing luminosity, and is probably 
isotropic and closely related to the accretion process of the central black hole for both radio- 
quiet and radio-loud sources. For radio-loud broad-line AGNs, the jet contribution to the 
soft X-ray emission is probably unimportant. The beaming effect due to the relativistic jet 
seems also weak in the soft X-ray band. Otherwise, we should expect the large scatter in 
the L x vs. L H p relation for radio-loud sources in our sample. 

The high redshift sources here usually have the measurement of the Mg II A2798 broad 
emission line instead of the H/? broad emission line in the observed wavelength. In Figure 9, 
we show the X-ray luminosity versus the Mg II emission line luminosity in our sample (left 
panel) and the radio luminosity versus the Mg II emission line luminosity (right panel). The 
result is consistent with what derived from Figure 8 for low redshift objects. 

3.4. Evolution of the black hole fundamental plane 

With the currently available large sample, we are able to investigate the evolution of the 
black hole fundamental plane relation by dividing the sample into different redshift bins. The 
result is listed in Table 4 and plotted in Figure 10. For radio-quiet objects, the coefficients 
^rm and ^r X are almost constant in different redshift bins. For radio-loud sources, ^rx is 
almost constant with redshift, while ^rm seems to decrease from positive to negative values 
as the redshift increases, with the exception of the last redshift bin, where £rm is positive 
again although it shows the largest error bars. 

Now we can compare the evolution result obtained by us with some theoretical predic- 
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tions. From Heinz & Sunyaev (2003), the radio flux produced via the synchrotron radiation 
from relativistic jet follows the scaling relation 

2p+13-(2+p)q r 2p+13+(p+6)ar 

F r oc M 8 +2p rh 2 (p+ 4 ) . (4) 

where a r and p represent the radio spectral index (/„ oc z/~ Qr ) and the electron energy 
distribution index (N oc E~ p ) respectively. M is the central black hole mass, rh is the 
dimensionless accretion rate (rh = - M ) . This relation is valid for radio-loud AGNs, whose 

V M E dd ' ' 

radio emission is believed to be produced via the synchrotron radiation of the relativistic jet. 

The dimensionless accretion rate (rh) is roughly proportional to the ratio between broad 
emission line and Eddington luminosities ( "^T^T^ ) (Wandel, Peterson & Malkan 1999). In 
the subsection above, we have shown that the soft X-ray luminosity is linearly scaled with the 
broad emission line luminosity for both radio-loud and radio-quiet subsamples. Therefore, 
we can use Lx/LEdd to represent the dimensionless accretion rate (rh). After replacing rh 
with L x /L Edd and considering the relation between the Eddington luminosity and the black 
hole mass (L Edd oc M BH ), equation (4) can be turned into 

2p+13+(p+6)a r 

L r (xM~ ar L x 2(p+4) . (5) 

Therefore, the coefficient £rm and £rx of the radio-loud subsample can be determined by 
a r and p. 

In our radio- loud AGN subsample, there are 114 sources that were also detected by the 
Green Bank 4.85GHz northern sky survey (GB6) (Gregore et al. 1996). Thus we are able 
to estimate the real radio spectral index a r (/„ oc z/~ Qr ) of this small subsample with the 
observed flux densities at 1.4GHz and 4.8GHz. We plot the radio spectral index versus the 
redshift for sources in this subsample in Figure 11. There is a weak trend that the spectral 
index increases with the increasing of the redshift. The mean values of a r in each redshift 
bins are -0.28, -0.12, -0.31, -0.27, -0.06, -0.11, and 0.04 (from low-z to high-z). We note 
that this is probably due to a selection effect in observations, as we may miss high redshift 
AGNs with the negative spectral index. 

Using the average radio spectral indeces in different redshift bins, we calculate the 
theoretical value of the coefficient £rm and £rx based on equation (5). The general trend of 
£rm is to decrease from 0.28 (a r = —0.28) to —0.04 (a r = 0.04), while £rx slightly increases 
from 1.23 to 1.44. This can be seen clearly from Figure 12. Therefore, the evolution of 
£rm and £ RX of the radio-loud subsample with the redshift is consistent with the theoretical 
prediction. 

On the other hand, it is also possible that the soft X-ray emission of broad-line AGNs is 
produced by the synchrotron process in a hot corona around accretion disks. Heinz (2004) 
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calculated the theoretical correlation between the radio luminosity, X-ray luminosity and 
black hole mass when the jet radiation dominates. If the physical parameters (i.e., magnetic 
field strength, electron energy distribution index, etc.) of the hot corona and the relativistic 
jet are similar, we can use the equations given by Heinz (2004) to roughly estimate the 
coefficients £rm and ^rx in different redshift bins. We find that with this model, 
decreases from 0.047 (a r = —0.28) to -0.007 (a r = 0.04), and £rx increases from 1.23 to 
1.44. Therefore, even if the soft X-ray emission mechanism is synchrotron radiation, the 
observed evolution of the black hole fundamental plane coefficients can still be explained. 
We can not exclude such a possibility because the origin of the soft X-ray emission is still 
uncertain. 

The radio spectral indeces of these 114 radio-loud AGNs obtained from observations 
allow us to directly calculate their rest-frame 1.4GHz radio fluxes without assuming the 
canonical value of a r as 0.5. The black hole fundamental plane relation is fitted based on this 
small subsample with measured a r . The derived £rx and £rm coefficients are 1.34 ±0.1 6 and 
—0.31 ± 0.21, respectively. These values are consistent within errors with what we obtained 
from the radio-loud subsample if we set a r = 0.5 (£rx = 1.50±0.08 and £rm = — 0.20±0.10). 
This indicates that our results for radio-loud AGNs are quite robust. 

4. DISCUSSION AND CONCLUSION 

4.1. Comparison with previous works 

The £rm coefficient we derived is rather small, meaning that the dependence of the 
fundamental plane relation on the black hole mass is weak, which is different from the result 
obtained by Merloni et al. (2003). However, there are some differences between our sample 
and that of Merloni et al. (2003). We only include the broad-line AGN, while their sample 
includes both GBHs and SMBHs. In our sample, the X-ray and radio emission are measured 
in 0.1-2.4keV and 1.4GHz respectively, while Merloni et al. (2003) used the data of 2-10keV 
X-ray emission and 5GHz radio core emission. It is still unclear whether the soft (0.1-2.4keV) 
and hard X-ray (2-10keV) emission have the similar origin for broad-line AGNs. A sample 
with both available data in the soft and hard X-ray bands will help us address this problem. 
This is beyond the scope of our current study and will be done in the near future. 

Another point worthy to mention is that the black hole mass data in Merloni et al. 
(2003) were mainly obtained from the literatures. This could introduce scatters due to the 
different mass measurement techniques. In our current study, the black hole masses of AGNs 
were estimated from the broad emission line and continuum properties by the same method, 



which does not introduce any additional bias. In addition, the black hole masses in our 
sample span a relatively smaller range than that in Merloni et al. (2003) since we only 
include broad-line AGNs. 



4.2. Relativistic beaming and other uncertainties 

Because the radio emission is produced by the relativistic jet, Doppler boosting of the 
synchrotron radiation (namely the relativistic beaming) would affect the observed radio 
flux significantly. Here we try to address the question that whether the observed larger 
fundamental plane coefficient £rx (or the larger radio luminosity) of radio-loud sources is 
due to the Doppler boosting effect. 

If the radio emission of both radio-loud and radio-quiet quasars are from jets, the larger 
radio luminosity observed in radio-loud sources can be considered to have much stronger 
Doppler boosting effect than radio-quiet AGNs. In other words, we can assume that the 
radio emission of radio-quiet AGNs is unbeamed. In subsection 3.3, we already show that 
the soft X-ray emission of broad-line AGNs is probably isotropic and unbeamed for both 
radio-loud and radio-quiet sources. Therefore, for a radio-loud quasar, its intrinsic radio 
luminosity (L r j et ) (unbeamed) may be estimated with its observed X-ray luminosity through 
the L r -L x correlation derived from radio-quiet sources. We will use the ratio between the 
observed radio luminosity (L r ) and the intrinsic radio luminosity (L r j et ) to represent the 
boosting factor of radio-loud sources. The equation of the Doppler boosting effect was given 
by Heinz k Merloni (2004): 



(1 + (3cos9) k+a - (1 - Pcos9) k+a * i' 

In Figure 13, we show the distribution of the boosting factor (L r /L r j et ) of the radio- 
loud subsample (left panel) and the boosting factor as a function of the inclination angle 
9 (right panel) when different Lorentz factor (r) is given. It is apparent in Figure 13 that 
only with smaller 9 (< 5°) and larger V (> 10), we can produce boosting factors as high as 
1000. Such conditions can only be met in BL Lac objects and are not likely to be the case 
for normal broad-line AGNs studied in this work. With the typical Lorentz factor of V ~ 5 
for broad-line AGNs (Orr & Browne 1982) and a non-negligible inclination angle (9 > 10°) 
(Maraschi et al. 1994), the boosting value is estimated to be less than 30. As can be seen in 
Figure 13, there are about half of radio- loud sources whose boosting factors are larger than 
the predicted maximum boosting value. Therefore, although the Doppler boosting effect 
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indeed has significant influence on the radio emission of radio-loud AGNs, this effect alone is 
not enough to explain the observed larger radio luminosity of radio-loud broad-line AGNs. 
This is consistent with the result given by Heinz & Merloni (2004) for the case of unbeamed 
X-ray emission. 

The non-simultaneous observations in the radio, X-ray and optical bands for sources 
in our sample may lead to other uncertainties. Both the ROSAT and FIRST surveys were 
conducted in the 1990s (Becker et al. 1995; Voges et al. 1999; Britzen et al. 2007). SDSS- 
I observations started from 1998, and ended in 2005 (York et al. 2000). So the data in 
different bands were obtained within ten years or so. Unless the luminosities of most objects 
in our sample varied significantly in these years, our result may not be affected too much 
due to this effect. In addition, several (2~3) factors change of the luminosities of some 
individual objects have little influence on the statistically significant results derived here for 
a large sample. However, as broad-line AGNs often show X-ray variabilities in timescales 
from hours to decades, the non- simultaneous observations could be an issue if the X-ray 
fluxes of AGNs vary significantly. Therefore, simultaneous observations in different bands 
for a sample of broad-line AGNs are still needed to confirm our results, although they are 
difficult to conduct for a large sample. 



4.3. Emission mechanisms 

The existence of the AGN radio-loud and radio-quiet dichotomy is still an unsolved 
issue. There are evidences that the radio emission of radio-quiet AGNs is likely produced 
by a weak (sub-relativistic) jet near the black hole (Blundell & Beasley 1998; Leipski et al. 
2006), while radio- loud quasars are usually associated with large scale jets of higher radio 
power (Rawlings & Saunders 1991; Miller, Rawlings & Saunders 1993). Here we fit the black 
hole fundamental plane relation for a broad-line AGN sample and find that the coefficients 
are quite different between radio-loud and radio-quiet quasars, especially the coefficient £,rx, 
which is steeper for radio-loud objects. We also find that the larger radio luminosities of 
radio-loud AGNs can not be produced by the Doppler boosting effect alone. Therefore, the 
radio emission mechanism may be quite different between radio-loud and radio-quiet AGNs. 
Such kind of difference can be caused by many uncertain physical parameters, such as the 
detailed disk/jet magnetic field strength and configuration, the electron energy distribution, 
and the black hole spin, etc. (Heinz & Sunyaev 2003; Sikora et al. 2007). More detailed 
observational and theoretical studies are still required. 

For the soft X-ray emission, we have shown that the correlation slopes between the soft 
X-ray luminosity and the broad emission line luminosity are all around 1 for both radio- 
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loud and radio-quiet AGNs. Therefore, the soft X-ray emission is probably produced via the 
accretion process near the central black hole, and it is unbeamed and isotropic for both radio- 
loud and radio-quiet sources. Otherwise larger scatters should exist in the Lx vs. Lh/j (L/Mgii) 
relation. However, there are many uncertain factors related to the origin of the soft X-ray 
emissions of AGNs (i.e., the accretion flow, hot corona, and warm absorber, etc.). A detailed 
study of the soft X-ray spectra for a larger sample of broad-line AGNs may give us clues to 
understand the origin of the soft X-ray emissions. 

4.4. Conclusion 

We revisited the fundamental plane relation of the black hole activity based on a large 
broad-line AGN sample selected on the basis of the cross-identification of the RASS, SDSS, 
and FIRST catalogs. The results of our work confirm the main result of Wang et al. (2006), 
namely, the black hole fundamental plane relation of the radio-quiet subsample is different 
from that of the radio-loud subsample; the coefficient £ RX becomes larger as the radio- 
loudness increases; the black hole mass seems unimportant in the black hole fundamental 
plane relation. We also found that the soft X-ray emission is most likely produced via 
the accretion process of the central black hole for both radio-quiet and radio-loud sources. 
In particular, for radio-loud sources, the jet contribution to the soft X-ray emission seems 
unimportant and the Doppler boosting of the relativistic jet is also weak in the soft X-ray 
band. Moreover, by dividing the radio-loud and radio-quiet samples into different redshift 
bins, we studied the evolution of the fundamental plane relation. For radio-quiet sources, 
there seems to be no clear evolution, while for radio-loud objects, the correlation coefficient 
£ RM tends to decrease as the redshift increases. We found that the evolution of the radio 
spectral index can help us at least partly understand such an evolution. Finally, we briefly 
discussed the beaming effect and some other uncertainties associated with the fundamental 
plane relation derived here. We found that Doppler boosting effect indeed has significant 
influence on the radio emission of radio-loud AGNs, but this effect alone is not enough to 
explain the observed radio luminosity of the radio-loud sources. 
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A. RADIO-LOUD NARROW-LINE SEYFERT 1 GALAXIES IN OUR 

SAMPLE 

Narrow-Line Seyfert 1 (NLS1) galaxies are a sub-class of AGN population. Their optical 
broad permitted emission lines are usually narrower (FWHMh/j < 2000 km/s ) than that in 
normal Broad-Line Seyfert 1 (BLS1) galaxies(Osterbrock & Pogge 1985). The NLSls also 
show weak [OIII]5007/H/5 to tai emission and strong Fell emissions (Boroson & Green 1992). 
Recently, Komossa et al.(2006) argued that the classical criteria (FWHMh/3 < 2000 km/s) 
is not well defined or even completely arbitrary. They suggested that -R4570 > 0.5 may be a 
physically more meaningful criteria to distinguish the NLS1 from the ordinary BLSlo 

The study of the physical mechanism of NLSls is still ongoing. There is growing evidence 
that most NLSls are objects with low black hole masses and high accretion rates, close to 
or even above the Eddington accretion rate (Collin & Kawaguchi 2004). 

Radio-loud NLSls are rare objects in NLS1 population (Komossa et al. 2006). In order 
to understand their radio properties, it is important to expand the number of radio-loud 
NLSls. Anderson et al. (2007) roughly examined the optical broad permitted emission line 
of AGNs in their catalog. They marked those objects that have FWHMh/3 < 2000 km/s as 
"NLS1?" in the comment columns of their tables. Seventy-four of these objects have optical 
spectra from SDSS and are detected in the FIRST radio survey, and thus are included in our 
analysis in this paper. We calculated the radio-loudness values for these objects, and found 
that five of them are radio-loud. We list the properties of the five sources in Table Al. 

In a newly published work, Yuan et al. (2008) present a comprehensive study of a 
sample of 23 radio-loud NLS1 galaxies. Among those 23 sources, two radio-loud NLSls 
are also discovered in our present work, which are SDSS J144318. 56+472556. 7 and SDSS 
J114654. 28+323652. 3. These two independent studies confirm the nature of radio-loud 
NLSls of these two objects. 



6 The optical Fe II strength, -R4570, is the ratio of the Fe II complex between the rest wavelength A4434A 
and A4684A to the total HP flux, including the narrow component (Boroson & Green 1992) 
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Table 1. The AGN Samples' 



Name z log(L x /l erg s" 1 ) log(L r /l erg s" 1 ) log R log(M/M ) L broad _ line b flag c 



SDSS J000608.04 - 010700.7 0.949 45.60 41.13 1.345 8.71 43.47 

SDSS J000710. 01 + 005329.1 0.316 44.95 39.79 0.549 9.07 42.90 1 

SDSS J004319.73 + 005115.4 0.308 44.64 39.61 0.506 9.42 42.90 1 



a Only 3 rows of the catalog are shown here. A complete catalog will be provided in the elecronic version 

b The broad emission line luminosity. It is represented by the Hp broad component luminosity or the Mg II emission line luminosity 
(when Hp is unavailable). 

c l means that Ljj^ is used for Lbroad— line: means that L/Mgii is used for Lt, r oad— ii nc . 



Table 2. Derived Fundamental Plane Relation 



Subsample 


Number 




£,RM 


Constant 




Total 


725 


1.47 ± 


0.06 


0.04 ±0.07 


-0.33 ± 0.06 


0.83 


Radio-quiet 


227 


0.73 ± 


0.10 


0.31 ± 0.12 


-0.68 ± 0.07 


0.42 


Radio-loud 


498 


1.50 ± 


0.08 


-0.20 ±0.10 


0.05 ± 0.10 


0.75 


Merloni et al. (2003) 




0.60 ± 


0.11 


n 7S +o.ii 


7 no+4.05 

'•°°-4.07 


0.88 


TotaP 


115 


1.33 ± 


0.15 


0.30 ±0.18 


-0.40 ±0.14 


0.89 


Radio-quict a 


39 


0.85 ± 


0.10 


0.12 ± 0.13 


-0.77 ±0.07 


0.38 


Radio-loud a 


76 


1.39 ± 


0.17 


0.17 ±0.21 


-0.17 ±0.21 


0.77 



a Results derived bv lWang et al . (2006) 



Table 3. Partial Correlation Test For The Fundamental Plane Correlation 



X 


Y 


Z 




Type 


Number 




T 












log Lx 


log Lr 


log 


D 


Radio-loud 


498 





321 





0264 


<1 


000E- 


10 


log Lx 


log Lr 


log 


D 


Radio-quiet 


227 





270 





0424 


1 


916E- 


10 


log Lx 


log Lr 


log 


D 


Total 


725 





308 





0220 


<1 


000E- 


10 


log Lx 


log Lr 


log 


M 


Radio-loud 


498 





567 





0201 


<1 


000E- 


10 


log Lx 


log Lr 


log 


M 


Radio-quiet 


227 





578 





0366 


<1 


000E- 


10 


log Lx 


log Lr 


log 


M 


Total 


725 





598 





0188 


<1 


000E- 


10 


log(Lx/Ledd) 


log(Lr/Ledd) 


log 


D 


Radio-loud 


498 





472 





0249 


<1 


000E- 


10 


log(Lx/Ledd) 


log(Lr/Ledd) 


log 


D 


Radio-quiet 


227 





517 





0340 


<1 


000E- 


10 


log(Lx/Ledd) 


log(Lr/Ledd) 


log 


D 


Total 


725 





412 





0197 


<1 


000E- 


10 


log(Lx/Ledd) 


log(Lr/Ledd) 


log 


M 


Radio-loud 


498 





572 





0204 


<1 


000E- 


10 


log(Lx/Ledd) 


log(Lr/Ledd) 


log 


M 


Radio-quiet 


227 





553 





0325 


<1 


000E- 


10 


log(Lx/Ledd) 


log(Lr/Ledd) 


log 


M 


Total 


725 





550 





0160 


<1 


000E- 


10 



Table 4. Derived Fundamental Plane Relation In Different Redshift Bins 



Rcdshift Range £rx £,rm N 



Radio-Quiet subsample 



2 < 0.13 


0.33 ±0.20 


0.23 ± 0.19 


77 


0.13 < z < 0.4 


0.35 ±0.23 


0.16 ± 0.20 


81 


0.4 < z < 1.0 


0.49 ±0.34 


0.12 ± 0.34 


48 


z > 1.0 


0.55 ±0.69 


0.45 ± 0.77 


21 


Radio- Loud subsample 


2 < 0.3 


0.63 ±0.24 


0.15 ± 0.23 


58 


0.3 < 2 < 0.45 


1.59 ±0.38 


-0.26 ± 0.24 


67 


0.45 < 2 < 0.65 


0.81 ±0.35 


-0.03 ± 0.21 


86 


0.65 < 2 < 0.80 


1.73 ±0.53 


-0.21 ± 0.24 


71 


0.80 < 2 < 1.0 


0.70 ±0.40 


-0.91 ± 0.41 


64 


1.0 < 2 < 1.3 


0.95 ±0.44 


-0.66 ± 0.34 


77 


2 > 1.3 


1.27 ±0.38 


0.11 ± 0.37 


75 



Table Al. Radio-Loud Narrow Line Seyfert 1 Objects 



Name 


z 


FWHM H(3 /1 km s" 1 


log(M BH /Mo) 


logR 


-R4570 


SDSS J144318.56 + 472556.7 a 


0.703 


1810.1 


7.14 


2.91 


5.50 


SDSS J114654.28 + 323652.3 


0.465 


2374.3 


7.43 


1.98 


1.59 


SDSS J073320.84 + 390505.2 


0.664 


2867.8 


7.86 


2.80 


3.27 


SDSS J154510.96 + 345246.9 


0.516 


3269.7 


8.13 


1.34 


0.74 


SDSS J123304.05 - 003134.1 a 


0.471 


3297.7 


7.81 


2.32 


2.18 



a _ff ( 3 emission line is fitted with just one Gaussian 
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Fig. 1. — Global properties of our AGN sample, with the top panels showing the histograms 
of the redshift and logarithm of radio-loudness, the middle panels showing the histograms of 
the X-ray luminosity and black hole mass (in units of M ), and the bottom panels showing 
the radio luminosity in logarithm units and ratios of X-ray to Eddington luminosity. 
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Fig. 2. — Rest-frame 1.4GHz radio luminosity (L r ) vs. the 0.1-2.4 keV X-ray luminosity, 
with different symbols representing different logarithmic bins of the black hole mass (in units 
of Mq). In the left panel we plot the logarithm of the luminosity directly, while in the right 
panel we scale the radio and X-ray luminosity with h-Edd- 
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Fig. 3. — Rest-frame 1.4GHz radio luminosity (L r ) vs. the 0.1-2.4 keV X-ray luminosity, 
with different symbols representing different radio- loudness bins. We plot the radio vs. X-ray 
luminosity in the left panel, and the Eddington luminosity scaled radio vs. X-ray luminosity 
in the right panel. 
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Fig. 4. — Ratio of the radio to X-ray luminosity vs. the black hole mass. Different symbols 
represent objects in different radio-loudness bins. 
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Fig. 5. — Edge-on view of the black hole fundamental plane for all AGNs in our sample. 
Different symbols represent sources in different radio-loudness bins. The dashed line is the 
best-fitting result for radio-quiet AGNs. 
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Fig. 6. — Correlation between log(L r / L Edd ) and 0.96 x log(L x /L Edd ). The symbols and 
lines have the same meanings as in Figure 5. 
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Fig. 7. — Results of the Monte Carlo test using scrambled radio fluxes. The position of the 
vertical line represents Pearson's correlation coefficient of the real dataset. 
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Fig. 8. — Plot of the 0.1-2.4keV X-ray luminosity versus the H/3 broad emission line lumi- 
nosity (left panel) and the rest-frame 1.4GHz radio luminosity versus the H/3 broad emission 
line luminosity (right panel). The meanings of symbols and lines are: green diamond: radio- 
loud AGNs; black cross: radio-quiet AGNs; black line: best-fitting result for the radio-quiet 
subsample; dashed black line: best-fitting result for the radio-loud subsample. 
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Fig. 9. — Plot of the 0.1-2.4keV X-ray luminosity versus the Mg II emission line luminosity 
(left panel) and the rest-frame 1.4GHz radio luminosity versus the Mg II emission line 
luminosity (right panel). The symbols and lines have the same meanings as in Figure 8. 
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Fig. 10. — Dependence of the black hole fundamental plane coefficients £rm and £rx on 
redshift for radio-quiet and radio-loud AGNs. 
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Fig. 11. — Redshift dependence of the radio spectral index a r of the 114 radio- loud AGNs 
detected both at 1.4GHz and 5GHz frequencies. 
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Fig. 12. — Predicted black hole fundamental plane coefficients based on equation (5). Here, 
a r is the radio spectral index and p equals to 2. 
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Fig. 13. — The histogram of the boosting factor (L r /L r j et ) of the radio-loud subsample (left 
panel) and the boosting factor as a function of the inclination angle 9 (right panel), with 
different line representing different Lorentz factor (r). 



